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^ '. Abstract. 
(N 

In this paper, we study negative order KdV (NKdV) equations and give their Hamiltonian struc- 

i tures. Lax pairs, infinitely many conservation laws, and explicit multi-soliton and multi-kink wave 

■ solutions thorough bilinear Backlund transformations. The NKdV equations studied in the paper 

are differential and can be derived from the first member in the negative order KdV hierarchy. The 

NKdV equations are not only gauge-equivalent to the Camassa-Holm equation through some hodo- 
(N ■ 

^ ' graph transformations, but also closely related to the Ermakov-Pinney systems and the Kupershmidt 

<N : 

I deformation. The bi-Hamiltonian structures and a Darboux transformation of the NKdV equations 

are constructed with the aid of trace identity and their Lax pairs, respectively. The 1- and 2- kink 
OO ' wave and soliton solutions are given in an explicit formula through the Darboux transformation. The 

o : . 

l-kink wave solution is expressed in the form of tanh while the 1-bell soliton is in the form of sech, 



and both forms are very standard. The collisions of 2-kink-wave and 2-bell-soliton solutions, are 
analyzed in details, and this singular interaction is a big difference from the regular KdV equation. 
Multi-dimensional binary Bell polynomials are employed to find bilinear formulation and Backlund 
transformations, which produce A^-soliton solutions. A direct and unifying scheme is proposed for ex- 
plicitly building up quasi-periodic wave solutions of the NKdV equations. Furthermore, the relations 
between quasi-periodic wave solutions and soliton solutions are clearly described. Finally, we show 
the quasi-periodic wave solution convergent to the soliton solution under some limit conditions. 
Keywords: Negative order KdV equations, bilinear Backlund transformation, Darboux 
transformation, kink wave solution, soliton solution, quasi-periodic solution. 
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10.5. Multi-periodic wave solutions 

1. Introduction 

The Korteweg-de Vries (KdV) equation 

Ut + Quu^ + u^^^ = 

was proposed by Korteweg and de Vries in fluid dynamics [38], starting from 
the observation and subsequent experiments by Russell [76]. There are many ex- 
cellent sources for the highly interesting background and historical development 
of the KdV equation, which brings it to the forefront of modern mathematical 
physics. In 1967, Gardner, Greener, Kruskal and Miura found the inverse scat- 
tering transformation method to solve the Cauchy problem of the KdV equation 
with sufficiently decaying initial data [18j. Soon thereafter. Lax explained the 
magical isospectral property of the time dependent family of Schrodinger opera- 
tors which is now called the Lax pair, and introduced the KdV hierarchy through 
a recursive procedure [13]. In the same year a sequence of infinitely many poly- 
nomial conservation laws were obtained with the help of Miura's transformation 

PUSH- 

There are tools to view the KdV equation as a completely integrable system 
by Gardner, and Zakharov and Faddeev [ST] . The bilinear derivative method 
was developed by Hirota to find A^-soliton solutions of the KdV equation [2B] . 
The KdV hierarchy was constructed by Lax [42j through a recursive approach, 
and further studied by Gel'fand and Dikii [2U]. On the base of the inverse spec- 
tral theory and algebro-geometric methods, the inverse scattering method was 
extended to periodic initial data by by Novikov, Dubrovin, Lax, Its, Matveev et 
al [TTl Hni EZ] • For more recent reviews on the KdV equation one may refer 
to literature [I1E1I1E1E1I2S1I1SIISS1IIZIEH]. 

All the work done in the above mentioned publications dealt with the positive 
order KdV hierarchy, which includes the KdV equation as a special member. 
However, there was only little work on the NKdV hierarchy. Verosky [80] studied 
symmetries and negative powers of recursion operator and gave the following 
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negative order KdV equation (called the NKdV equation thereafter) 

Vt = w^, 

(1.1) 

Wxxx + 4:VWx + 2VxW = 0. 

and Lou [47J presented additional symmetries based on the invertible recursion 
operator of the KdV system and particularly provided the following NKdV equa- 
tion (called the NKdV-1 equation thereafter) 

vt = 2uux, + f M = 0, ( — ) +2uux = 0, (1.2) 

\ u /t 

which can be reduced from the NKdV equation (1.1) under the following trans- 
formation 

w = u , V = . (1.3) 

u 

Moreover, the second part of NKdV-1 equation (1.2) is a linear Schrodinger equa- 
tion or Hill equation 

Uxx + VU = 0. 

Fuchssteiner [17] pointed out the gauge-equivalent relation between the NKdV 
equation (1.1) and the Camassa-Holm (CH) equation [6] 

i^t + fTLxU + 2mUx = 0, m = u — u^x 

through some hodograph transformation, and later on Hone proposed the as- 
sociate CH equation, which is actually equivalent to the NKdV equation (1.1), 
and gave soliton solutions through the KdV system [52] • Zhou [85j generalized 
the Kupershmidt deformation and proposed a kind of the mixed KdV hierarchy, 
which contains the NKdV equation (1.1) as a special case. 

Very recently, Qiao and Li [5T] gave a unifying formulation of the Lax repre- 
sentations for both negative and positive order KdV hierarchies, and furthermore 
studied all possible traveling wave solutions, including soliton, kink wave, and pe- 
riodic wave solutions, of the integrable NKdV-1 equation (1.2) with the following 
Lax pair 

Lip = ipxx + vip = Xip, 

(1.4) 

tpt = ^M^A - huu^X V- 
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The most interesting [61] is: the NKdV-1 equation has both sohton and kink 
solutions, which is the first integrable example, within our knowledge, having 
such a property in soliton theory. 

Studying negative order integrable hierarchies plays an important role in the 
theory of peaked soliton (peakon) and cusp soliton (cuspon). For instance, the 
well-known CH peakon equation is actually produced through its negative order 
hierarchy while its positive order hierarchy includes the remarkable Harry-Dym 
type equation [02] ■ The Degasperis-Procesi (DP) peakon equation [9] can also 
be generated through its negative order hierarchy [M]. Both the CH equation 
and the DP equation are typical integrable peakon and cupson systems with 
nonlinear quadratic terms [SJ dUl SSI E21 ES]- Recently, some nonlinear cubic 
integrable equations have also been found to have peakon and cupson solutions 

[321 ESI ESI ES]. 

In this paper, we study the NKdV hierarchy, in particular, focus on the 
NKdV equation (1.1) and the NKdV-1 equation (1.2). Actually, as per [iTflSS]. 
the NKdV equation (1.1) can embrace other possible differential- integro forms 
according to the kernel of operator K = jd^ + ^{vdx + dxv). Here we just list 
the NKdV-1 equation (1.2) as it is differential and also equivalent to a nonlinear 
quartic integrable system: 

UUxxt — UxxUt — 2u^Ux = 0. 

The purpose of this paper is to investigate integrable properties, A^-soliton 
and A^-kink solutions of the NKdV equation (1.1) and NKdV-1 equation (1.2). In 
section 2, the trace identity technique is employed to construct the bi-Hamiltonian 
structures of the NKdV hierarchy. In section 3, we show that the NKdV equa- 
tion (1.1) is related to the Kupershmidt deformation and the Ermakov-Pinney 
systems, and is also able to reduced to the NKdV-1 equation (1.2) under a trans- 
formation. In section 4, a Darboux transformation of the NKdV equation (1.1) 
is provided with the help of its Lax pairs. In section 5, as a direct application of 
the Darboux transformation, the kink-wave and bell soliton solutions are explic- 
itly given, and the collision of two soliton solutions is analyzed in detail through 
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two-solitons. In section 6, an extra auxiliary variable is introduced to bilinearize 
the NKdV equation (1.1) through binary Bell polynomials. In section 7, the bi- 
linear Backlund transformations are obtained and Lax pairs are also recovered. 
In section 8, we will give a kind of Darboux covariant Lax pair, and in section 
9, infinitely many conservation laws of the NKdV equation (1.1) are presented 
through its Lax equation and a generalized Miura transformation. All conserved 
densities and fluxes are recursively given in an explicit formula. In sections 10, 
a direct and unifying scheme is proposed for building up quasi-periodic wave 
solutions of the NKdV equation (1.1) in an explicit formula. Furthermore, the 
relations between quasi-periodic wave solutions and soliton solutions are clearly 
described. Finally, we show the quasi-periodic wave solution convergent to the 
soliton solution under some limit conditions. 

2. Hamiltonian structures of the NKdV hierarchy 

To find the Hamiltonian structures of the NKdV hierarchy, let us re-derive 
the NKdV hierarchy in matrix form. 

2.1. The NKdV hierarchy 

Consider the Schrodinger-KdV spectral problem 



where A is an eigenvalue, ijj is the eigenfunction corresponding to the eigenvalue 
A, and ^; is a potential function. 

Let (fii — ip, (fi2 — i^x-i then the spectral problem (2.1) becomes 



where ip = ((^1,^92)'^ is a two-dimensional vector of eigenfunctions. 

The Gateaux derivative of spectral operator U in direction ^ at point v is 



(2.1) 




(2.2) 




(2.3) 



which is injective and linear with respect to the variable ^. 
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The Lenard recursive sequence {Gm} of the spectral problem (2.1) is defined 



by 

G^i e KerK = {G\KG = 0}, Go E KerJ = {G\JG = 0} 

KGm-i = JG^m, m = 0,-l,-2--- , 
which directly produces the NKdV hierarchy: 



(2.4) 



vt = KGm-i = JGm, m = -l,-2--- (2.5) 

where 

K = ^dl + ^{vd^ + d^,v), J = d^, (2.6) 
and K is exactly a recursion operator of the well-known KdV hierarchy 

Vt = K''v^, n = 0,l,2,---. 

The first equation (m = 0) in the NKdV hierarchy (2.5) is trivial equation 

vt = JGo = 0, JGo = KG-i = 0. 

The second equation (m = —1) in the NKdV hierarchy (2.5) takes 

vt = KG-i = 0, 

which is exactly the NKdV equation (1.1) by replacing G^i = w. 

In a similar way to the paper [61], we construct zero curvature representation 
for NKdV hierarchy. 

Proposition 1. Let U be the spectral matrix defined in (2.2), then for an 
arbitrarily smooth function G G C°°(]R), the following operator equation 

- [U, V] = U'[KG] - XU'lJG] (2.7) 

admits a matrix solution 



V = ViG) 



( -'-G. '-g\ 

4 2 

1 111 

\ Grr VG H AG —Gr I 

\ 4 2 2 4 / 



which is a linear function with respect to G, and Gateaux derivative is defined 
by (2.3). 

Theorem 1. Suppose that {Gj, j = — 1,— 2,---} is the first Lenard se- 
quence defined by (2.4), and Vj = V{Gj) is a corresponding solution to the 



operator equation (2.7) for G — Gj. With Vj being its coefficients, a mth matrix 
polynomial in A is constructed as follows 

m 

Then we conclude that the NKdV hierarchy (2.5) admits zero curvature repre- 
sentation 

Ut-Wm,, + [U,W^] = 0, 

which is equivalent to 

(2-8) 

i=i \-^Gj^xx - 2'^Gj + ^XGj jGj^xJ 
This theorem actually provides an unified formula of the Lax pairs for the whole 
NKdV hierarchy (2.5). 

According to theorem 1, the NKdV equation (1.1) admits Lax pair with 
parameter A 

Lip = ipxx + vip^ Xip, 

or equivalently, 

= {dl + v)il} = A^, 

(2.9) 

Mi) = {Adldt + Avdt + 2wd^ + ?>w^)ip = 0. 
The NKdV equation (1.1) also possesses Lax pair without parameter 

2.10 

MtP = {Adldt + Avdt + 2wd^ + 2>w^)ij = 0. 

Especially, taking the constraint v — —Uxx/u and w — v? & KerK, we then 
further get the NKdV equation (1.2) and its Lax parir (1.4). 

2.2. Hamiltonian structures 



Proposition 2. [77] For the spectral problem (2.2), assume that is a 
solution to the following stationary zero curvature equation with the given ho- 
mogeneous rank 

V,= [U,V]. (2.11) 
Then there exists a constant /3, such that 

holds, where (■, ■) stands for the trace of the product of two matrices. 

Let {Cm, = —I, — 2 ■ ■ ■ } be the negative order Lenard sequence recursively 
given through (2.4) and 

-1 



m=—oo 



be a series with respect to A. Assume that V\ = V{Gx) is the matrix solution for 

the operator equation (2.9) corresponding to G = Gx. So, V\ can be written as 

-1 

m=—oo 

Then, we have the following proposition. 
Proposition 3. Vx satisfies the following Lax form 

Vx,. = [U,Vx]. 

Proof. By (2.4), we have 

-1 -1 
{K-XJ)Gx= Y KG^^-^- Y ^G'„A-™+i 

m=— oo m=—OD 
-1 

= KG^.X-' + Y iKG„,-i - JGm)X-^ = 0. 

m=— oo 

Therefore, Proposition 1 implies 

Vx,. - [U, Vx] = U'[KGx\ - XU'[JGx\ = U'[KGx - XJGx] = 0. 

□ 

Next, we discuss the Hamiltonian structures of the hierarchy (2.5). It is 
crucial to find infinitely many conserved densities. 
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Theorem 2. 

(1) The hierarchy (2.5) possesses the bi-Hamiltonian structures 

Vt^K—^—^J-^, m = -l,-2---, (2.14) 

where the Hamiltonian functions are imphcitly given through the 
following formulas 

H_i^G^ie KerK, //^ = — , m = -l,-2---. (2.15) 

m 

(2) The hierarchy (2.5) is intcgrable in the Liouville sense. 

(3) The Hamiltonian functions are conserved densities of the whole 
hierarchy (2.5) and therefore they are in involution in pairs. 

Proof. A direction calculation leads to 

By using the trace identity (2.12) and the expansion (2.13), we obtain 
J2 Gm\-A= J] (m-l-/3)G'^_iA— + (-l-/3)G'_i, 



S_ 

Sv 



\m=—oo 



(2.16) 



m = -l,-2--- . 

If taking G_i 7^ 0, form (2.16) we find /3 — —1 and 

5Hm 



--Gm-1: m = -l,-2---, (2.17) 

where if^ are given by (2.15). Substituting (2.17) into (2.5) yields the bi- 
Hamiltonian structures (2.14). 

Next, we consider infinitely many conserved densities to guarantee integra- 
bihty of the hierarchy (2.16). Since J and K are skew-symmetric operators, we 
infer that 

C*J = {J-^Kyj - -K* =JC, 

which implies 



6v 6v 

= {C^+'G^i, JC^-^Go) = {i^n+i, Hm-i}, m,n< -1. 
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Repeating the above argument gives 

{Hn,Hm} = {Hm,Hn} = {Hm+n, H^i}. (2.18) 

On the other hand, we find 

{H„,, Hr,} = (C^G^,, J^G^r) = {rC^'G^,, C'G^,) = HJ^. (2.19) 

Then combining (2.18) with (2.19) leads to 

{Hm, Hn} = 0, 

which imphes that {Hm} are in involution, and therefore the hierarchy (2.14) are 
integrable in Liouville sense. 

Especially, under the constraint (1.3), we obtain bi-Hamilton structures of 
the NKdV equation (1.2) 

vt = K— — = J^, 
ov ou 

where two Hamiltonian functions are given by 

Ho = ^u^, H_i = -u^, 

which can also be written in a conserved density form in the sense of equivalence 
class 

Ho ~ — - J u^dx, H^i ^ ~ J "^^dx. 

3. Relations to other important equations 
3.1. Kupershmidt deformation 

Recently a class of new integrable systems, known as the Kupershmidt de- 
formation of soliton equations, have attracted much attention. This topic starts 
from Kupershmidt, Karasu-Kalkani' work [271 [STf [39]. 

For the Lenard operator pair (2.6), we define Lenard gradients recursively by 

KGj = JGj+i, KG^i = JGo = 0, j = 0, ±1, ±2, ■ ■ • , 
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then KdV hierarchy is 

Vt = KGm-i = JGm, m = 0,±l,±2,--- (3.1) 

which contains both the NKdV hierarchy and the positive order KdV hierarchy. 
The first equation (m = 0) in the KdV hierarchy (3.1) is trivial system 

= JGo = 0, KG^i = JGo = 0, (3.2) 

which can be regarded "sharp threshold" equation of the NKdV hierarchy 

and positive order KdV hierarchy. 

A Kupershmidt nonholonomic deformation of the hierarchy (3.1) takes 
Vt = JGm + Jw, m = 0, ±1, ±2, ■ ■ ■ , 

(3.3) 

Kw = 0, 

where two operators K and J are given by (1.4). Then the first flow (m = 0) of 
the hierarchy (3.3) is exactly the NKdV equation (1.1) 

Vt = w^, 

WxXX + 4:VWx + 2VxW = 0, 

which may be regarded as a Kupershmidt nonholonomic deformation of the 
threshold equation (3.2) 

3.2. NKdV hierarchy with self-consistent sources 



Soliton equations with self-consistent sources have important physical appli- 
cations, for example, the KdV equation with self-consistent source describes the 
interaction of long and short capillary-gravity waves [H| \5U\ EH E2] . 

For the N distinct Xj of the spectral problem (2.1), the functional gradient 
of Xj with respect to v is 



6v 

Here we define the whole KdV hierarchy with self-consistent sources as follows 

6X ^ 

Vt = JGm + aJ-r- = JGm + ClJ/^lph 

OV ^-^ ■' 

(3.4) 

1pj,xx + {V + Xj)ll)j = 0, 

m = 0,±l,±2,--- ; j = l,---,iV. 
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Taking m = 1 in the hierarchy (3.4) leads to the KdV equation with self- 
consistent sources 

N 



1 ^ 

= j{vxxx + Qvvx) + adx ^ iJj, 



4 

'>Pj,xx + {v + Xj)ipj = 0, i = 1, • • • , A^, 
while choosing m = — 1 in the hierarchy (3.4) gives the NKdV equation with 
self-consistent sources 

N 

Wxxx + ^VWx + 2VxW = 0, 

'il^j,xx + {v + Xj)'ilJj = 0, j = 1, • • • , A^. 
Obviously, taking N — 1, m — 0, a — 1, v ^ v + Xi in the hierarchy (3.4), 
then we get the NKdV equation (1.2) 

Vt = ii^Dx, 'ipl,xx + Vi^i = 0, 

which may be regarded as the threshold equation (3.2) with self-consistent sources. 
3.3. Reduction of the NKdV equation (1.1) 

Theorem 3. (u, v) is a solution of NKdV-1 equation (1.2) if and only if (ty, v) 
with w — is a solution of NKdV equation (1.1) under the transformation 

Uxx + vu = 0, (3.5) 

which is actually a linear Schrodinger equation or Hill equation. 
Proof. Let 

w = u^, (3.6) 

then by (1.1), we have 

Vt^Wx^ '^UUx, 

which is the first equation of (1.2). By (3.6), the second equation of (1.1) leads 
to 

^Ux{Uxx + VU) + u{Uxx + Vu)x ^0, 
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or equivalently, 

[u^{u^x + vu)]cc = 0, (3.7) 

Apparently, according to (3.7), if {u,v) is a solution of the NKdV-1 equation 
(1.2), then {w,v) is a solution of the NKdV equation (1.1) where w = v?. Re- 
versely, if {w^v) is a solution of the NKdV equation (1.1), then {u,v) is also a 
solution of the NKdV-1 equation (1.2) under the transformation (3.5). 

For a given function 0, let us define the following Baker- Akhiezer function 

u = exp (pdcc^ , (3.8) 
then (3.2) yields the following Riccati equation 

(j)x + (f)'^ + v = 0. (3.9) 

So, we have 

Theorem 4. {u,v) is a solution of the NKdV-1 equation (1.2) if and only 
if {w, v) is a solution of the NKdV equation (1.1) as is a solution of the Riccati 
equation (3.9) while u is the Baker-Akhiezer function (3.8) and w = u^. 

3.4. Ermakov-Pinney equation 



The Ermakov-Ray-Reid systems 

1 



1 '0 
(l>xx + W^(x)0 = —-^G{ — ), 



were originally introduced by Ermakov [70] . Due to their nice mathematical 
properties of Ermakov systems admitting a novel integral of motion together 
with a concomitant nonlinear superposition principle and extensively physical 
applications, there has been numerous an extensive literature devoting to the 
analysis of the Ermakov systems EOl EH [Z21 E31 El] • The most simple case is 
equation 
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which is called the Ermakov-Pinney equation. The Ermakov-Pinney equation 
is a quite famous example of a nonlinear ordinary differential equation. Such 
an equation (and generalizations thereof) have been shown to be relevant to a 
number of physical contexts including quantum cosmology, quantum field theory, 
nonlinear elasticity and nonlinear optics [l6l [751 EH] ■ A recent account of some 
of its properties along with applications in cosmological settings can be found in 

Ref. [sg. 

Proposition 4. Suppose that {w,v) is a solution of the NKdV equation 
(1.1). Let 

w = pt = ij'^, V =px, 

then ip satisfies a Ermakov-Pinney equation 

ij.. + v^ = ^, (3.9) 

where fi is an integration constant. 

Especially, if {u,v) is the solution of the NKdV-1 equation (1.2), let 



u = 0exp yi J jj,(f)~ dxj , (3.10) 
then satisfies the Ermakov-Pinney equation 

+ W0 = (3.11) 

Proof. Substituting transformation w = ip'^ into the second equation of the 
NKdV equation (1.1) yields 

Wxxx + ^VWx + 2VxW = 2ilj{ip^^ + V1p)x + Qlpxilpxx + vip) 

2 

= -^[{i^xx + vij)ij% = 0, 

which leads to (3.9). 

Substituting transformation (3.10) into the second equation of the NKdV 
equation (1.2) yields 

Ux + vu = {(pxx + vcp — ^) exp j ii(j)~^dx^ = 0, 

which implies (3.11). 
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Proposition 5. For given function v, let ■0i)'02 are two solutions of linear 
Schrodinger equation 

Uxx + vu^O, (3.12) 

then equation 

Wxxx + 4:WxV + 2vxW = (3.13) 

admits a general solution 

w^aipl + 26ViV'2 + cipl (3.14) 

where 

OC - 6^ = W = 1pl1p2,x - '01,x'02- 

Proof. Let w — ip'^, hy proposition 4, then ip satisfies a Ermakov-Pinney 
equation (3.9). It is easy to check that if ipi and ■02 are two solutions of equation 
(3.14), then 

ip = ^ aipf + 2bipiip2 + cip2 
is a solution of equation (3.9). So (3.14) is a general solution of the equation 
(3.13). 

4. Darboux transformation of NKdV equations 

In this section, we shall construct a Darboux transformation for general 
NKdV equation (1.1), and then reduce it to the NKdV-1 equation (1.2). 

4.1. Dcirboux transformation 

A Darboux transformation is actually a special gauge transformation 

Tp^Tijj (4.1) 

of solutions of the Lax pair (2.9), here T is a differential operator (For the Lax 
pair (2.10), the Darboux transformation with A = can be obtained). It requires 
that ijj also satisfies the same Lax pair (2.9) with some L and M, i. e. 

Ltp = A^, L = TLT-\ 

(4.2) 
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Apparently, we have 

which imphes that L and M are required to have the same forms as L and M, 
respectively, in order to make system (2.9) invariant under the gauge transforma- 
tion (3.4). At the same time the old potentials u and f in L, M will be mapped 
into new potentials u and v in L, M. This process can be done continually and 
usually it may yield a series of multi-soliton solutions. 

Let us now set up a Darboux transformation for the system (2.9). Let ■0o = 
ipo{x,t) be a basic solution of Lax pair (2.9) for Aq, and use it to define the 
following gauge transformation 

= T^, (4.3) 

where 

T^d^-a, a = d^lni/;o- (4.4) 
Prom (2.9) and (4.4), one can see that a satisfies 

aa: + a'^ + v-X = (4.5) 

4:crxxt + i'^cr-^crt + Avut + 2w(7^ + QaUxt + 3^^^ = 0. (4.6) 
Proposition 6. The operator L determined by (4.2) has the same form as 
L, that is, 

L = dl + v, 

where the transformation between v and v is given by 

{i = V + 2(7j;. (4.7) 

The transformation: (^,f) — )■ (^,w) is called a Darboux transformation of the 
first spectral problem of Lax pair (2.9). 

Proof. According to (4.2), we just prove 

LT = TL, 

that is, 

{dl + v){d^-a) = {dx-a){dl + v), 
which is true through (4.5) and (4.7). 
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Proposition 7. Under the transformation (4.3), the operator M determined 
by (4.2) has the same form as M, that is, 

M = Adldt + Avdt - 2wd^ - 3w^, (4.8) 

where the transformations between w, v and w, v are given by 

w = w + 2(Tt, V = V + 2ax. (4.9) 

The transformation: {ip,w,v) {ip,w,v) is Darboux transformation of the sec- 
ond spectral problem of Lax pair (2.9). 

Proof. To see that M has the form (4.8) same as M, we just prove 

MT = TM, (4.10) 

where 

M ^ Adldt + fdt + gd, + h, (4.11) 

with three functions f,g, and h to be determined. Substituting M, M, L into 
(4.10) and comparing the coefficients of all distinct operators lead to: 
coefficient of operator dxdt 

/ = 4v + 8(7a; = Av, 

which holds by using (4.9). 
coefficient of operator 

g ^2w + Aat^ 2w, 

which implies from (4.9). 
coefficient of operator dx 

h = Sa^t + 5wx — 2aw + ga = 6axt + Sw^; + 2{ax + cr^ + v)t 
= QOxt + "iwx = 3Wx, 
here we have used equation (4.5) and (4.9). 
coefficient of operator dt 

—A(7xx ~ — 4t'a; — Ava, 

that is, 

+ 2aax + ■Ua; = 0. 
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which holds by using (4.5). 
coefficient of non-operator: 

4:crxxt + f<^t + g<^x + crh + 3w^^ - 3aw^ = 0, 

that is, 

4a^^t + Ua^at + ivat + 2wa^ + Gaa^t + 3w^^ = 0, 

which is the equation (4.6). We complete the proof. □ 

Propositions 4 and 5 tell us that the transformations (4.3) and (4.9) send the 
Lax pair (2.9) to another Lax pair (4.2) in the same type. Therefore, both of the 
Lax pairs lead to the same NKdV equation (1.1). So, we call the transformation 
{ip,w,v) — >■ {ip,w,v) a Darboux transformation of the NKdV equation (1.1). In 
summary, we arrive at the following theorem. 

Theorem 5. A solution w, v of the NKdV equation (1.1) is mapped into 
its new solution w, v under the Darboux transformations (4.3) and (4.9). 

4.2. Reduction of Dcirboux transformations 

To get Darboux transformations for the NKdV-1 equation (1.2), we consider 
two reductions of the Darboux transformations (4.3) and (4.9). 

Corollary 1. Let X = > 0, then under the constraints w = u'^, v = 
—Uxx/u, the Darboux transformations (4.3) and (4.9) are reduced to a Darboux 
transformation of the NKdV-1 equation (1.2) {il!,v,u) {ip,v,u), where 

^i^^T^jj, 'D = i; + 2(7^, u^k-^{u:,-au) ^k-^Tu. (4.13) 

Proof. For A > 0, suppose that {v,u) is a solution of the NKdV-1 equation 
and -0 is an eigenfunction of the Lax pair (1.4), then we have 

Therefore, the Lax pair (1.4) can be rewritten as 

fpxx + v'iIj^ Xi/j, 

1 _i 1 _i (4-12) 

ipt = -uX {wpx - u^ij) = -ud^ {wp) = N{u, X)i), 

where N — N{u, X) — \ud~^u. 
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According to Proposition 6, the spectral problem of Lax pair (4.12) is covari- 
ant under the transformation (4.13), that is, 

So, we only need to prove 

7Pt^N{u,X)7P (4.14) 
Substituting (4.13) into the left hand side of (4.14) yields 

= ^[{'^x - (Tu)d-\uilj) - ^(7V(^ix - au)d-\u'ilJo)] (4.15) 
— -ku[d~^ {uijj) + k~'^{ux — au)il)]. 

In the same way, substituting (4.13) into the right hand side of (4.14) produces 

N{u, = ^ud~^[k~^{u^ - au){2p^ - aip)] 

= ^k~'^u[u^'ilj - d~^{u^^'ilj) - aui/j + 9j^(V'(rVo,a;xiiV')] (4-16) 
— -A;~"^'u[A;^9~"^(mV') + {ux — (Tu)iI)]. 

Combining (4.15) with (4.16) imphes that (4.14) holds. □ 
In a similar way, we also have the following result. 

Corollary 2. Let A = 0, then under the constraints w = u^, v = —Uxx/u, 
the Darboux transformations (4.3) and (4.9) are reduced to another Darboux 
transformation of the NKdV-1 equation (1.2) {ip,v,u) — )> {■ip,v,u), where 

v^v + 2ax, -0 = V - '0o'^<^^^^(V'o'0), 
u — 'iljQ^ad~^{'ijjou), uj^O. 



u — 



with (7 = a^ln(l + a-lV'o)■ 



5. Applications of the Darboux transformation 
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In this section, we shall apply the Darboux transformations (4.3) and (4.9) to 
obtain kink- type and bell- type of exphcit solutions for the NKdV equation (1.1). 



5.1. The kink- wave solutions 

For the case of A = A;^ > 0, we substitute v — 0,w — 1 into the Lax pair (2.9) 
and choose the following basic solution 

ijj = + e~^ = 2 cosh^, ^ = kx-^t + 'y, (5.1) 

where 7 and k are two arbitrary constants. 
Taking kf, then (4.4) and (5.1) lead to 

(7i a^^lnV^ = /citanh^i, ^1 = kix - -^t + ^i. 

zki 

The Darboux transformation (4.9) gives bell-type solution for the NKdV equation 
(1.1) 

— 2(7i X — 2A;?sech^^i, 

w = 1 — 2ai^t = tanh 4i- 
By using Darboux trasformation (4.13), we get a kink- type wave solution for the 
NKdV equation (1.2) 

^ k^^{ux - au) ^ -tanh ^i, ^ - + -fi. (5.3) 

Remark 1. There is much difference between traveling waves of the NKdV 
equation (1.2) and of the classical KdV equation. For the NKdV equation (1.2), 
its one-wave solution is a negative-moving (i.e. from right to left) kink-wave with 
velocity — 1/2/c^, amplitude ±1 and width 1/ki. Its amplitude is independent of 
velocity, and width is directly proportional to the velocity. For the KdV equation 

ut + 6uux + Uxxx = 0, (5.4) 

one-soliton solution is 

2 2' ^ ' 

which is a bell-type positive-moving wave with velocity /c^, amplitude k'^/2 and 

width 1/k, respectively. Its amplitude is directly proportional to velocity, and 

width is inversely proportional to the velocity. 
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(a) 



(b) 





Figure 1. The two-kink wave solution u{x,t) with parameters: ki = 1, 
k2 = 0.6. (a) Perspective view of the wave, (b) Overhead view of the wave, 
with contour plot shown. The bright lines are crests and the dark lines are 
troughs. 

Let US now construct two-kink solutions to see the interaction of two kink 
solutions. According to (4.4), 

^ = T^/' = (a,-ai)(e« + e-«) 



is also an eigenfunction of Lax pair (2.9). Taking A = we have 



0"2 = —ki tanh^i 



_ ^2 

"-1 "-2 



(5.6) 



(5.7) 



ki tanh — k2 tanh ^2 
Repeating the Darboux transformation (4.9) one more time, we get two soli- 
ton solution for the the NKdV equation (1.1) 

[kf - fc|)(fc|sech2^2 - /t^sech^^i) 



v^^ = v^ + 2a2,x 



2a. 



{ki tanh ^1 — k2 tanh ^2 
ki tanh ^2 — ^2 tanh ^1 



2d 



ki tanh — k2 tanh ^2 , 
Therefore, we obtain a two-kink wave solution of the NKdV equation (1.2) 

~ A;2 tanh ^1 - A;i tanh ^2 „x 

u = — ; — -. (5.8) 

Kitanh^i — A;2tanh^2 

Let us use the two-kink wave solution (5.8) to analyze interaction of the two 
one-soliton solutions. Without loss of generality, we suppose ki > k2 > 0, then 
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(a) 



(b) 



(c) 



(d) 



(e) 




Figure 2. interaction between singular soliton csch^i and smooth soliton 
sech^2 with parameters: (a) t = —3, (b) t = —0.05, (c) t = 0, [d) t — 0.05, (e) 
t ^ 3. 



we have 



Therefore, on the fixed fine =constant, we get 



tanh ^2 



-1, t — )■ +00, 



and it follows (5.8) that 

k2 tanh^i + ki 



u ~ — = coth \ t,i In — 

ki tanh^i + k2 \ 2 ki + k2 



t +00. 



(5.9) 



In a similar way, one can get 

tanh ^2 ~ 1 as t — > — oo, 
which are main parts compared with terms 1 and e^^\ and it follows (3.19) that 

, t ^ -oo. (5.10) 



u 



In a similar way, on the line ^2 =constant, we will arrive at 

, as t ^ +00, 



u ~ tanh I ^2 + - In ^ ^ 



2 ki + k2 



^~tanh(6-^ln|-p^,, 



as t ^ —00. 



(5.11) 
(5.12) 
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Remcirk 2. Prom expressions (5.9)-(5.12), we see that the two-kink wave 
solution (5.8) is a singular solution, which is able to be decomposed into a singular 
kink-type solution and a smooth kink wave solutions. The expressions (5.10) and 
(5.12) show that the wave tanh^2 is on the left of the wave coth^i before their 
interaction, while the expressions (5.9) and (5.11) show that the wave coth^i is 
on the left of the wave tanh ,^2 after their interaction. The shapes of the two kink 
waves coth ^1 and tanh ^2 don't change except their phases. Their phases of the 
two waves coth^i and tanh {2 are In > and — In < 0, respectively as 
the wave is negatively going along the x— axis. Very interesting case is particular 
ai t — 0: collision of such two kink waves forms a smooth bell-type soliton and 
its singularity disappears (See Figure 2). 

After their interaction It can be seen that the the two kink waves resume 
their original shapes. At the right moment of interaction, the two kink waves are 
fused into a smooth bell-type soliton. The two-kink wave interactions possess the 
regular elastic-collision features and pass through each other, and their shapes 
keep unchanged with a phase shift after the interaction. Here, wc also demon- 
strate a fact that the large-amplitude kink wave with faster velocity overtakes 
the small-amplitude one, after collision, the smaller one is left behind. 

5.2. The bell-type soliton solutions 

(i) For the case of A = (i.e. without parameter A), we substitute v = 
—k^, w — into the Lax pair (2.10), and choose the following basic solution as 

ip = + e~^, ^ = kx + -^t, 

2k 

where k is an arbitrary constant. 
Taking k = ki, (4.4) gives 

a = ai = dxlnip = /citanh^i, ^1 = kix + -^rj—t. (5.13) 
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Using the Darboux transformation (4.9), we have one-sohton solution for the 
NKdV equation (1.1) 

V = V + 2o"i X = 2/c^sech^^i — k^, 

(5.14) 

w — — 2(Ji,t = sech^^i. 

So, we get a one-sohton solution for the NKdV-1 equation (1.2) by using 
Darboux transformation (4.17) 

u — sech^i, ^1 = kix + -^t. (5.15) 

Remark 3. For the negative order KdV equation (1.2), its one-soiton solu- 
tion (5.15) is a smooth bell-type negative-moving wave, whose velocity, amplitude 
and width are l/2kl, ±1 and 1/ki, respectively. Its amplitude is independent of 
velocity, and width is directly proportional to the velocity. 

(ii) For the case of A = — A;^, we take a seed solution oi v — — 2A;^, ty = 1 in 
the Lax pair (2.9), and choose the following basic solution as 

where k is an arbitrary constant. 
Taking k = ki sends (3.7) to 

(7 = (7i = In-j/) = /ci tanh^i, ii = kix — ——t + ^i. (5.13) 

Using the Darboux transformation (3.12), we then get one-soliton solution 

= v + 2ai^j; = -2A;itanh^^i + 71, 

(5.14) 

= l-2ait^ 1 + sech^^i. 



which cannot satisfies the constraint (3.3), so vw^ is not soliton for the NKdV 
equation (1.2). 

Remark 4. For the NKdV equation (1.1), its one-soiton solution (5.14) is a 
smooth bell-type positive-moving wave, whose velocity, amplitude and width are 
l/2kl, ±1 and l/ki, respectively. Its amplitude is independent of velocity, and 
width is directly proportional to the velocity. 
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Let's construct a two-soliton solution of the NKdV equation (1.1). According 
to the gauge transformation (4.4), 

^ = T^ = {d,-a,)ie^ + e-^) 

is also an eigenfunction of Lax (2.9). We have 

kf — ki 

^ ^ ~'~ /ci tanh^i — A;2 tanh^2 

Repeating the Darboux transformation (4.9) one more time, we obtain 

-// -/ , o (fc? - fci)(A;isech^^2 - fcjsech^^i) 

' (fc,tanhei-fc2tanh6)^ ' 

_ 2^ _ / /ci tanh ^2 - ^2 tanh ^1 V 
\A;i tanh^i — A;2 tanh^2/ ' 
which is the same for NKdV-1 equation (1.2). So we get two-soliton solution with 

(5.8) 

~ ^ki tanh ^2 — ^2 tanh ^1 
ki tanh^i — k2 tanh ^2 ' 
but here = kjX — ^t, j = 1, 2. 

6. Bilinearization of the NKdV equation 

The bilinear derivative method, developed by Hirota [2E] , has become a pow- 
erful approach to construct exact solutions of nonlinear equations. Once a non- 
linear equation is written in a bilinear form by using some transformation, then 
multi-solitary wave solutions or quasi-periodic wave solutions can usually be ob- 
tained [23 EDI EH ESI SSI EH [13] • However, unfortunately, this method is not as 
direct as many people might wish because the original equation is reduced to two 
or more bilinear equations under new variables called Hirota's variables. Since no 
a general rule to select Hirota's variables, there is no rule to choose some essen- 
tial formulas (such as exchange formulas), either. Especially the construction of 
bilinear Backlund transformation relies on a particular skill and appropriate ex- 
change formulas. On the other hand, in recent years Lambert and his co-workers 
have found a kind of the generalized Bell polynomials playing important role in 
seeking the characterization of bilinearized equations. Based on the Bell polyno- 
mials, they presented an alternative procedure to obtain parameter families of a 
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bilinear Backlund transformation and Lax pairs for soliton equations in a quick 
and short way SHI SI]. 

6.1. Multi-dimensional binary Bell polynomials 

The main tool we use here is a class of generalized multi-dimensional binary 
Bell polynomials [26]-|41]. 

Definition 1. Let > 0, k = I,-- - ,i denote arbitrary integers, / = 
f{xi, ■ ■ ■ be a C°° multi-variable function, then 



=exp(-/)9,"^-- 9,"; exp(/) (6.1) 

is a polynomial in the partial derivatives of / with respect to a;i, ■ ■ ■ ,Xi, which 
we call a multi-dimensional Bell polynomial (a generalized Bell polynomial or 
F-polynomial). 

For the two dimensional case, let / = f{x,t), then the associated Bell poly- 
nomials through (6.1) can produce the following representatives: 

^xif) = fx, ^2z(/) = f2x + fx, ^3x(/) = fsx + 3/x/2x + fx, 
^x',i(/) = fx,t + fxft, Y2x,t{f) = f2x,t + f2xft + '^fx,tfx + fx ft, ' ' ' ■ 

Definition 2. Based on the use of above Bell polynomials (6.21), the 
multi-dimensional binary Bell polynomials ( 3^-polynomials) are defined by 



'ynixi,--- jHixti^d , ^) ^nixi,--- ,nixi(^f^ 



ri + ■ ■ ■ + re is odd. 



.hr^xu-,rexi, \-re is even, 

which is a multi-variable polynomial with respect to all partial derivatives g-nxi,--- ,rixi 
(ri H h r£ is odd) and /irixi,-,r^x^ (^i H h is even), = 0, ■ ■ ■ , n^, k = 

o,---,i. 

The binary Bell polynomials also inherit partial structures of the Bell poly- 
nomials. The first few lower order binary Bell Polynomials are 

yxig) = Qx, y2x{g, h) = h2x + gl, yxAa, h) = hxt + gxgt- 

(6.2) 

y^xig, h) = g3x + 3gxh2x + gl,--- ■ 
Proposition 8. The link between binary Bell polynomials y-mxi,--- ,nexi{g, h) 
and the standard Hirota bilinear expression D^^ ■ ■ ■ D^^F ■ G can be given by an 



(6.5) 
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identity 

yn,.,,-,n,.A9 = Ini^/C, h = InFG) = {FCy'D^l ■ ■ ■ Dl^F ■ G, (6.3) 

whereni + n2 + ■ • • + > 1, and operators D^^^, • • • , D^^ are classical Hirota's 
bilinear operators defined by 

D^l ■ --D^jF-G = (a,,-a,0"^ • ■ ■ (9.,-9..)"^F(a:i, • • • ,xe)G{x[, ■ ■ ■ ,4)1.,=,,,..,,,=,,. 

In the special case of F = G, the formula (6.4) becomes 

F-'D2l ■ ■ ■ D^G ■ G = :^n,.,,..,„,.,(0, ? = 21nG) 

0, rii H h is odd, (6.4) 

Pn ixi,- ,nixt{Q)-i Til -\- ■ ■ ■ -\- n(, is even. 
The first few P-polynomial are 

P2x{q) = q2x, PxM) = (ixu P^xiq) = q^x + '^qlx, 
P&x{q) = 96a; + '^^q2xqAx + is^L, • • • . 

The formulas (6.4) and (6.5) will prove particularly useful in connecting nonlinear 
equations to their corresponding bilinear forms. This means that if a nonlinear 
equation is expressedby a linear combination of P-polynomials, then the nonlinear 
equation can be transformed into a linear equation. 

Proposition 9. The binary Bell polynomials ymxi,- ,rnxi{v-iw) can be sepa- 
rated into P-polynomials and F-polynomials 

{FG)~'^Dll ■ • • D^jF ■ G = ynixu-,nexe{9, h)\g=inF/G,h=lnFG 

= ynixi,-,nexe{9,9 + q,)\g=lnF/G,q=2lnG ^, 

(6.6) 

n-iH \-n(=even ri=0 r(=0 i=l ^ ^ 

The key property of the multi-dimensional Bell polynomials 

implies that the binary Bell polynomials y-nixi,--- ,nixe{g, h) can still be linearized 
by means of the Hopf-Cole transformation g = Inip, that is, ip = F/G. The 
formulas (6.6) and (6.7) will then provide the shortest way to the associated Lax 
system of nonlinear equations. 
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6.2. Bilinecirization 

Theorem 6. Under the transformation 

V = Vo + 2{hiG)2x, w = Wo + 2{hiG)xt, 

the NKdV equation (1.1) can be bihnearized into 

(L*^ + UvoDl - D,Dy)G • G = 0, 

(6.8) 

{2DtDl + GwoDl + DtDy)G ■ G = 0. 
where y ia an auxihary variable, and Uq, Vq are two constant solutions of the 
NKdV equation (1.1). 

Proof. The invariance of the NKdV equation (1.1) under the scale transfor- 
mation 

X ^ Xx, A"i, V X~'^v, w X~°'~^w 

shows that the dimensions of the fields v and w are —2 and — (« + 1), respectively. 
So we may introduce a dimensionless potential field q by setting 

v^Vo + q2x, w^Wq- q^f (6.9) 

Substituting the transformation (6.9) into the equation (1.1), we can write the 
resulting equation in the following form 

q4x,t + M2xq2x,t + 2?3xfe + ^Voq2x,t + 2Wofe = 0, 

which is regrouped as follows 
2 1 

g?4x,t + 2(fe?2x,t + QxtqSx) + :^q4x,t + 2?2x?2a;,t + ^Voq2x,t + '^WoqSx = 0, (6.10) 

where we will see that Such an expression is necessary to get a bihnear form of the 
equation (1.1). Further integrating the equation (6.10) with respect to x yields 

2 1 
E{q) = ^{q3x,t + ^2xqxt + Swofe) + :^d~^dt{q4x + ^qlx + 12^^og2x) = 0. 

(6.11) 

In order to write the equation (6.11) in a local bilinear form, let us first get 
rid of the integral operator d~^. To do so, we introduce an auxiliary variable y 
and impose a subsidiary constraint condition 

q4x + + 12^^og2a; " qxy = 0. (6.12) 
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Then, the equation (6.10) becomes 

H(l3x,t + ^q2xqxt + 3Wog2x) + Qyt = 0. (6.13) 

According to the formula (6.5), the equations (6.12) and (6.13) are then cast 
into a pair of equations in the form of P-polynomials 

PUq) + - Pxy{q) = 0, 

^PsxAq) + ^WoP2x{q) + Pyt{q) + 37 = 0. 
Finally, by the property (6.4), making the following variable 

q^2hiG <^ ^; = 't;o + 2(lnG')2a;, f/; = + 2(ln G')^;*, 

change above system to the following bilinear forms of the NKdV equation (1.1) 
as follows 

{Dl + l2v^Dl D,Dy)G -0 = 0, 

(6.14) 

{2DtDl + GwoDl + DtDy)G ■ G = 0, 
which is also simultaneously bilinear system in y. This system is easily solved 
with multi-soliton solutions by using the Hirota's bilinear method. □ 

Finally, we show that the NKdV-1 equation (1.1) can be directly bihnearized 
through a transformation, not Bell polynomials. Making dependent variable 
transformation 

V = vo + 2{lnF%^, u = G/F, (6.15) 
we can change the equation (1.2) into 

2{F^t - F^Ft) = G\ 

FxxG — 2FxGx + GxxF + vqFG = 0, 
which is equivalent to the bilinear form 

D,DtF-F = G\ {Dl + vo)F-G = 0. (6.16) 

It is obvious that the bilinear form of the NKdV-1 (6.16) is more simple than the 
bilinear form of NKdV (6.15). 

6.3. N-soliton solutions 
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As usual as the normal perturbation method, let us expand G in the power 
series of a small parameter e as follows 

G'=l + £^(l)+£V')+£3^(3) + ... 

Substituting the above equation into (6.7) and sorting each order of £, we have 
£ : {Dl + UvoDl - D,Dy)g, -1 = 0, 

(6.16) 

{2DtDl + GwoDl + DtDy)g^'^ -1 = 0, 
: {Dt + 12voDl - D,Dy){2g^^^ ■ 1 + g^"^ ■ g^'^) = 0, 

(6.17) 

(2DtDl + QwoDl + DtDy){2g^^^ ■ 1 + g^^^ ■ g^^^) = 0, 

: {Dt + UvoDl - D^Dy){g^'^ • 1 + ^« • g^'^) = 0, 

{2DtDl + GwoDl + DtDy){g^^^ ■ 1 + ^^^^ • g^^^) = 0, (6.18) 



By employing formulae mentioned above, the system (6.16) is equivalent to 
the following linear system 



which has solution 



^a)=e«, ^^kx-^^^t + ik' + 12vok)y + a, (6.19) 

where k and a are two arbitrary parameters. 

Substituting (6.12) into (6.10) and (6.11) and choosing g^"^") — g^^^ — ■•■ — (), 
then the G's expansion is truncated with a finite sum as 

G = 1 + 

which gives regular one-soliton solution of the NKdV equation (1.1) 
v^Vo + 2dl ln(l + e^) = + ysech2C/2, 

w = Wo + 2dtd, ln(l + e«) = + — ^sech^^/S, (6.20) 
where 7 = (/c^ + 12vok)y + a, and k, Vq, Wq are constants. 
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Let Wo — 1, vo — 0, then the solution (6.20) reads as a kink- type solution of 
the NKdV-I equation (1.2) 

2 

ii = ±tanh^/2, ^ — kx — —t + ^. 
In a similar way, taking 

^(i) = e«^+e^., i^ = kjX-^^^^t + ^„ J = 1,2, 

we get a two-soliton wave solution 

v = vo + 2dl ln(l + + e^^ + e^^+^^+^i^) 



t« = Wo - 29^4 ln(l + + + e^^+^^+^i^)^ 

/ — ACo \ ^ 

>li2 = In 



/ci — k 



(6.21) 



. ^1 + ^2 , 

In general, we can get a N-soliton solution of the NKdV equation (1.1) 

(N N 

^lj=Q,l j=i i<j<N 

(N N 

Mj=o,i j=i i<i<iv 

2 



An = In 



/v/ 

where the notation X]^j=o i represents all possible combinations = 0, 1, and 

6- = /^V^-|S^ + ^^-' J = l,2,---,iV. 

In the following, we discuss the soliton solutions for NKdV-1 equation by 
using bilinear equation (6.16). Let us expand F and G in the power series of a 
small parameter e as follows 

(2)^2 , .(4)^4 , .(6)^6 



F=l + f'>e' + + f^^S^ + 



Substituting the above equation into (6.16) and arranging each order of £, we 
have 

^i'i + ^0^^'^ = 0, 

6.23 



2/i? = (9^'r, 

(6.24 



Let Vq — —k"^, it follows from the first equation of (6.23) and (6.24) that 

Substituting (6.25) into the second equation of (6.23) leads to 

- k^9^'^ = 0, 

from which we may take g^'^^ — 0, further choose g^^^ = • • • = 0, f^^^ = • • • = 0. 
So F and G are truncated with a finite sum as 

4 

Finally, the formula (6.14) gives one-soliton solution of the NKdV-1 equation 
(1.2) 

V — 2/c^sech^^ — k"^, u — sech^. 



7. Bilinecir Backlund transformation 

In this section, we search for the bilinear Backlund transformation and Lax 
pair of the NKdV equation (1.1). 

7.1. Bilinecir Backlund transformation 
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Theorem 7. Suppose that F is a solution of the bilinear equation (6.8), and 
if G satisfying 

(L>^-A)F.G = 0, 

(7.1) 

[DtDl + 2woD,, + (4t;o + 3A) A]F -(3 = 0, 
then G is another solution of the equation (6.8). 
Proof. Let 

be two different solutions of the equation (6.10). Introducing two new variables 

h^{q + q)/2 = In(FG), g ^ {q - q)/2 = ln(F/G'), 

makes the function E invariant under the two fields q and q: 
E{q)-E{q)^E{h + g)-E{h-g) 

= ^vog^t + 4:Wog2x + '^Qzx^t + 4:h2xgx,t + 4:K,t92x + 4:d~^{h2xg2x,t + h2x,t92x) 
= 2d,{y2x,t{9, h) + ^vM9) - 2wo3^.(5)) + R{9. h) = 0, 

(7.2) 

where 

R{g, h) = -2d^[{h2x + 9l)9t] + ^h2xgxt - ^h2x,t9x + 49^^(^2xfl'2x,t + h2x,t92x)- 

This two-field invariant condition can be regarded as a natural ansatz for a bilinear 
Backlund transformation and may produce some required transformations under 
additional appropriate constraints. 

In order to decouple the two-field condition (7.2), let us impose a constraint 
so as to express R{g, h) in the form of x-derivative of 3^-polynomials. The simple 
choice of the constraint may be 

y2x{9. h) = h2x + 9l-\ (7.3) 

which directly leads to 

R{g, h) = 2Xg^t + ^h2xgxt - 4:h2x,t9x - ^9l9xt = ^>^9xt, (7.4) 
where h2x,t — ~'^9x9xt and h2x — ^~ 9x used. 
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Using the relations (7.2)-(7.4), we derived a coupled system of iV-polynomials 

(7.5) 

y2.,t{9, h) + (4t;o + ^\)yt{g) + 2w^y^M = 0, 
where we prefer the second equation to be expressed in the form of conserved 
quantity without integration with respect to x. This is very useful to construct 
conservation laws. Apparently, the identity (6.2) directly sends the system (7.5) 
to the following bilinear Backlund transformation 
(L>^-A)F-G = 0, 

(7.6) 

[DtDl + 2woD^ + (4t;o + 3A) A]F ■ G = 0, 
where we have integrated the second equation in the system (7.5) with respect 
to X, and wq is the corresponding integration constant. □ 

7.2. Inverse scattering formulation 

Theorem 8. The NKdV equation (1.1) admits a Lax pair 

(7.7) 

Proof. By the transformation v = ln'0, it follows from the formulas (6.5) and 
(6.6) that 

ytig) = V't/V', yx{g) = V'x/V', y2x{g: h) = q2x + V'2a;/V', 

y2x,t{9^ h) = ^q^ti'x/i^ + q2xipt/ip + i^2x,t/i^, 

which make the system (7.5) linearized into a Lax pair with parameter A 

LV'=(9^ + g2.)V' = AV', (7.8) 

M^|; = [dtdl + {Avo + q2x)dt + 2(fe + wo)d, + SA^^jV^, (7.9) 
or equivalently, 

i^2x +V1p^ Xtp, 

4V'2x,t + 4f V't - ^wipx - ^w^iIj = 0, 
where the equation (7.8) is used to get the second equation. One can easily verify 
from equations (7.8) and (7.9) that 

[L, M] = g4^,t + 4(^0 + q2x)q2x,t + 2g3x(fe + Wo) = 
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exactly gives the NKdV equation (1.1) through replacing vq + q2x and wq + Qxt 
by V and w, respectively. □ 

8. Dcirboux covariant Letx pair 

In this section, we will give a kind of Darboux covariant Lax pair, whose form 
is invariant under the gauge transformation (4.3). 

Theorem 9. The NKdV equation (1.1) possesses the following Darboux 
covariant Lax pair 

Lip = Xip, 

McovV'-O, M^oy = M + 3dxL, 
under the gauge transformation tjj — Ti/j. This is actually equivalent to the Lax 
pair (2.9). 

Proof. In section 4, we have shown that the gauge transformation (4.1) maps 
the operator L{q) onto a similar operator 

L{q) = TL{q)T-\ 

which satisfies the following covariance condition 

L{q) — L{q + Aq) , q — q + Aq, with Ag = 21n0. 

Next, we want to find a third order operator Mco^iq) with appropriate coeffi- 
cients, such that Mcoviq) is mapped by gauge transformation (8.1) onto a similar 
operator Mcov{q), which satisfies the covariance condition 

Mcov(g) = M,ov(g + Aq), q^q + Aq. 

Suppose that is a solution of the following Lax pair 

Lip = Xip, 

(8.2) 

McovV' = 0, Mcov = ^dtdl + hdx + bidt + 63, 
where 61, &2 and 63 are functions to be determined. Then, the transformation T 
is required to map the operator M^ov to the similar one 

TM,^^T-^ = Meov, ^2,cov = 4=dtdl + hd^ + 62^* + h, (8.3) 
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where bi, 62 and 63 satisfy the covariant condition 

b, = bj{q) + Abj = bj{q + Aq), j = 1, 2, 3. (8.4) 
It follows from (8.2) and (5.3) that 

Abi = h-b^=4at, A62 = 62 - &2 = 8a,, (8.5) 

A63 = 63 - 63 = aAbi + 8a^t + bi,^,, (8.6) 

and a satisfy 

4(T2a;,t + ha^ + kat + aAbs + 63,, = 0. (8.7) 

According to the relation (8.4), it remains to determine 61, 62 and 63 in the form 
of polynomial expressions in terms of g's derivatives 

= 9a;, Qy, Qxy, (l2x, (l2y, q2x,y, ' ' ■), j = 1, 2, 3 

such that 

AFj = Fj{q + Aq, q^ + Ag,, qt + Aqt,---)- Fj{q, q^, qt,---) = Abj, (8.8) 

with Aqkx,it — 2(ln0)jtj;^/t, k,l — 1, 2, • • • , and Abj being given through the rela- 
tions (8.5)-(8.7). 

Expanding the left hand of the equation (8.8), we obtain 

Abi = AFi = F.^gAq + F.^g^Aq, + F.^^^Aqy + F,,,^^Aq,t + ■ ■ ■ = 4at = 2Aq,t, 

which implies that we can determine bi up to a arbitrary constant ci, namely, 

bi = Fi{qxt) = 2qxt + Ci, Ci is an arbitrary constant (8.9) 

Proceeding in the same way deduce the function 62 as follows 

b2 = F2{q2a^) = iq2:c + C2, (8.10) 

where C2 is an arbitrary constant. 

We see from the relation (8.6) that A63 contains the term bi^x — Q2x,t, which 
should be eliminated such that A63 admits the form (8.8). By the Lax pair (8.2), 
we have the following relation 

q2x,t = -CTxt - 2aat. (8.11) 
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Substituting (8.9) and (8.11) into (8.6) yields 

A63 = Aaat + 8a^t + '^q2x,t = Q(^xt = 3Ag2x,t- 

If choosing 

b3 = F3{q2^,t) = ^2x,t + C3, (8.12) 

the third condition 

AF3 = F3,,Ag + Fa,,^ Ag^ + Fa^.^Ag^ • • • = A63 

holds, where C3 is an arbitrary constant. 

Letting Ci = —2vq, C2 = 0, C3 = Wq in (8.9), (8.10) and (8.12), then it follows 
from (8.2) that we have the following Darboux covariant evolution equation 

McoviIj = 0, Mcov = 4:dtdl + 2q^tdx + 4:q2xdt + 3q2x,t: 

which coincides with the equation (8.7). Moreover, the relation between two 
operators L2,cov and L2 are related through 

The compatibility condition of the Darboux covariant Lax pair (8.2) exactly gives 
the NKdV equation(l.l) in Lax representation 

[Mcov, L] = g4^,t + 4(^0 + q2x)(l2x,t + (l3x{.(lxt + Wo) 
= Vxxx + '^VWx + 2VxW = 0. 

□ 

In the above repeated procedure, we are able to obtain higher order opera- 
tors, which are also Darboux covariant with respect to T, to produce higher order 
members of the negative order KdV hierarchy. 

9. Conservation laws of NKdV equations 

In this section, we will present infinitely many conservation laws in a local 
form for the NKdV equation (1.1) based on a generalized Miura transformation. 
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Theorem 10. The NKdV equation (1.1) possesses the following infinitely 
many conservation laws 

i=^n,t + G,,. = 0, n = l,2,-- - . (9.1) 

where the conversed densities F^s are recursively given by recursion formulas 
explicitly 

_ _^ ^ _ (9.2) 

fe=0 k=0 

and the fluxes G'^^s are 

Go = 2w/o = 2wv, Gi = 2wli — —2wVx, 

(9.3) 

Gn = 2wJ„, n = 2, 3, • • • . 
Proof. For the simplicity, let us select = ■u;o = in the transformation 
(6.9). We introduce a new potential function 

q2x^V + ^Vx + £W, (9.4) 
where e is a constant parameter. Substituting (9.4) into the Lax equation (7.10) 
leads to 

= [L, M] = (1 + ed^ + 2e^r))[-A{r) + eW)r]t - 2(?x - eri)tr]^ + r)2x,tl 

which implies that v = q2x, w = q^t given by (9.4) are a solution of the NKdV 
equation (1.1) if satisfies the following equation 

-4(77 + e^rf)r]t - 2{q^ - er))tr)x + V2x,t - ^Vt = 0. (9.5) 
On the other hand, it follows from (9.5) that 

[{Qx - er])t\x ^ -{V + eW)t- 
Therefore, the equation (9.5) can be rewritten as 

{ri2x - rf)t + [2r]{e^r] - q:c)t]x = 0, 
or a divergent-type form 

{r)2. + 2e^r]r], - r)% + {2r)w), = (9.6) 

by replacing q^t^w- 
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Inserting the expansion 

oo 

77 = •••)£", (9-7) 

n=0 

into the equation (9.4) and comparing the coefficients for power of £, we obtain 
the recursion relations to calculate in an explicit form 

-^0 = q2x ^v, -Io,x = -Vx, 

_ _ (9.8) 

— ~In-l,x ~ IkIn-2-ki H — 2, 3, • • • . 

fc=0 

Substituting (9.7) into (9.6) and simplifying terms in the power of e provide 
us infinitely many conservation laws 

Fn,t + Gn,x = 0, n = 1, 2, • • • 

where the conversed densities F^s and the fluxes G^s are by (9.2) and (9.3), 
respectively. □ 

Here, we already give recursion formulas (9.7) and (9.8) to show how to gener- 
ate conservation laws (9.6) based on the first few explicitly provided. Apparently, 
the first equation in conservation laws (9.6) 

'^xxt ~ 2l>l't + 2wVx + 'iWxV — 

is exactly the NKdV equation (1.1) 

Vt + Wx^ 0, 

'^XXX + ^VWrc + 2WVx — 0. 

which is reduced to the NKdV equation (1.2) under the constraints v = —Uxx/u 
and w = V?. 

In conclusion, the NKdV equation (1.1) is completely integrable and admits 
bilinear Backlund transformation. Lax pair and infinitely many local conservation 
laws. 



10. Quasi-periodic solutions of the NKdV equation 
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In this section, we study quasi-periodic wave solutions of the NKdV equa- 
tion (1.1) by using bihnear Backlund transformation (7.1) and bihnear formulas 
derived in section 9. 

In fact, a quasi-periodic solution, also called algebro-geometric solutions or fi- 
nite gap solutions, was originally studied in the late 1970s by Novikov, Dubrovin, 
McKean, Lax, Its, and Matveev et al [HI [361 SSI E^, based on the inverse 
spectral theory and algebro-geometric method. In recent years, this theory has 
been extended to a large class of nonlinear integrable equations including sine- 
Gordon equation, Camassa-Holm equation, Thirring model equation, Kadomtsev- 
Petviashvili equation, Ablowitz-Ladik lattice, and Toda lattice [3 ED [221 [221 [2S1 
I21[2I1[SZ1I521[SS1IS21[HZ1ISS]. The algebro-geometric theory, however, needs 
Lax pairs and is also involved in complicated analysis procedure on the Riemann 
surfaces. It is rather difficult to directly determine the characteristic parameters 
of waves, such as frequencies and phase shifts for a function with given wave- 
numbers and amplitudes. On the other hand, the bilinear derivative method 
developed by Hirota is a powerful approach for constructing exact solution of 
nonlinear equations in an explicit form. If a nonlinear equation is able to be writ- 
ten in a bilinear form by a dependent variable transformation, then multi-solitary 
wave solutions are usually obtained for the equation (211 [301 1311 ISSl [81] . Based on 
the Hirota forms, Nakamura proposed a convenient way to find a kind of explicit 
quasi-periodic solutions of nonlinear equations [56], where the periodic wave so- 
lutions of the KdV equation and the Boussinesq equation were obtained. Such 
a method indeed displays some advantages over algebro-geometric methods. For 
example, it does not need any Lax pair and Riemann surface for the given non- 
linear equation, and is also able to find the explicit construction of multi-periodic 
wave solutions. The method relies on the existence of the Hirota's bilinear form 
as well as arbitrary parameters appearing in Riemann matrix [131 HI] • 

10.1. Multi-dimensional Riemann theta functions 



42 

Let us first begin with some preliminary work about multi- dimensional Rie- 
mann theta functions and their quasi-periodicity. The multi- dimensional Rie- 
mann theta function is defined by 

t9(C,s, s|r)= ^ exp{27rz(C + £,n + s) - 7r(T(n + s),n + s)}, (10.1) 

where n = (rii, ■ ■ ■ , um)'^ G is an integer value vector, and s = (si, ■ ■ ■ , sn)'^, s = 
(ei, ■ ■ ■ ,eiy)'^ G is a complex parameter vector. <^ = (Ci, ■ " " ^Cn)'^, Cj = 
ajX + f3jt + 6j, aj,f3j,6j E Aq, j = 1,2, ■■■ ,iV are complex phase variables, 
where x, t are ordinary real variables and 6' is a Grassmann variable. The inner 
product of two vectors / = (/i, ■ ■ ■ , /at)"^ and g = {gi, ■ ■ ■ , Qn)'^ is defined by 

{f,9) = fiQi + 1292 H h fNQN- 

The matrix r = (rjj) is a positive definite and real- valued symmetric N x N ma- 
trix. The entries Tij of the periodic matrix r can be considered as free parameters 
of the theta function (10.1). 

In this paper, we choose r to be purely imaginary matrix to make the theta 
function (10.1) real-valued. In definition (10.1) for the case of s = e = 0, we 
denote i){^,T) = 'i9(<^, 0, 0|t) for simplicity. Therefore, we have 'i9(<^, s, 0|t) = 
^(C + s,t). 

Remark 4. The above periodic matrix r is different from the one in the 
algebro-geometric approach discussed in [57j-[llj, where it is usually constructed 
on a compact Riemann surface F with genus N eN. One may see that the entries 
in the matrix r are not free and difficult to be explicitly given. □ 

Definition 3. A function g{x,t) on x C is said to be quasi-periodic in 
t with fundamental periods Ti, ■ ■ ■ , G C if Ti, ■ ■ ■ , are linearly dependent 
over Z and there exists a function G{x,y) G x C'^ such that 

Gix, 7/1, ■ ■ ■ , Vj+Tj, ■■■ ,yk) = G{x, yi, ■ ■ ■ ,yj, ■ ■ ■ , yk), for all G C, j = 1, ■ ■ ■ , k. 

G{x,t,--- ,t, ■■■ ,t) = g{x,t). 

In particular, g{x,t) becomes periodic with T if and only if Tj = rrijT. □ 
Let's first see periodicity of the theta function 'i9(<^, r). 
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Proposition 10. [55] Let ej be the j— th column of N x N identity matrix 
In] Tj be the j— th column of r, and Tjj the (j,j)-entry of r. Then the theta 
function r) has the periodic properties 

■&{C + + ^Tj' ''") = exp(-27riCj + TTTjj)'&{C, ''')■ 

The theta function d{C,, r) which satisfies the condition (5.4) is called a mul- 
tiplicative function. We regard the vectors {e^, j = 1, ■ ■ ■ , A^} and {iTj, j = 
1, ■ ■ ■ , A^} as periods of the theta function {}{(^, r) with multipliers 1 and exp{—2'Ki(j + nrjj), 
respectively. Here, only the first vectors are actually periods of the theta func- 
tion t), but the last A^ vectors are the periods of the functions (9|^ ^ In r) 
and ln[^9(C + e, t)/^{C + h,T)], k,l = l,--- ,N. 

Proposition 11. Let Cj and Tj be defined as above proposition 2. The 
meromorphic functions /(<^) follow 

{^^) /(0 = %ln^|^|^, C, e, h G C^ j = l,---,Ar. 
then in all two cases (i) and (ii), it holds that 

f{C + ej+tTj) = f{C), CeC^^ j = l,---,A^, 

which implies that /(<^) is a quasi-periodic function. 

10.2. Bilinear formulae of theta functions 

To construct a kind of explicitly quasi-periodic solutions of the NKdV equa- 
tion (1.1), we propose some important bilinear formulas of multi-dimensional 
Riemann theta functions, whose derivations are similar to the case of super bi- 
linear equations [15j, so we just list them without proofs. 

Theorem 11. Suppose that '(9(<^, e', 0|t) and '(?(<^, e, 0|t) are two Riemann 
theta functions, in which £ = {ei, . . . ^En), s' = (e'^, . . . , e'^y), and C, = (Ci, ■ ■ ■ , Ca^); 
(j = OjX + Ujt + Sj, j = 1,2, ■ ■ ■ , N. Then operators Dr^, Df and 5* exhibit the 
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following perfect properties when they act on a pair of theta functions 
D,i9iC,e',0\T)-i»iC,e,0\T) 

= d.m, ^' - ^, -M/2|2r) |c=o^(2C, e' + £, ix/2\2t), (^O-^) 

where = (yUi, ■ ■ ■ ,/iAr), and the notation represents 2^ different transfor- 
mations corresponding to all possible combinations = 0, 1; ■ ■ ■ ; = 0, 1. 

In general, for a polynomial operator H^D^, Dt) with respect to and Dt, 
we have the following useful formula 

H{D,, A)^(C, 0|t) • ^(C, £, 0|t) = J] C(e', e, /u)^(2C, e' + e, /x/2|2t), 

(10.3) 

in which, explicitly 

C(e, e', /u) = ^ exp [-27r(T(n - /u/2), n - /Lt/2) - 27ri(n - /u/2, e' - e)] 

(10.4) 

where we denote A4 = {Am{n — ^/2, a), ATii{n — /i/2, a;)). 

Remcirk 6. The formulae (10.3) and (10.4) show that if the following equa- 
tions are satisfied 

C(e, e', /u) = 0, (10.5) 

for all possible combinations /Xi = 0, 1; //2 = 0, 1; • ■ ■ ; jijq = 0, 1, in other word, 
all such combinations are solutions of equation (10.5), then ■d{<^,e',0\T) and 
■j?(C, e, 0|t) are A^-periodic wave solutions of the bilinear equation 

H{D,, DtMC, s', 0|r) • i?(C, £, 0|t) = 0. 

We call the formula (10.5) constraint equations, whose number is 2-^. This for- 
mula actually provides us an unified approach to construct multi-periodic wave 
solutions for supersymmetric equations. Once a super symmetric equation is writ- 
ten bilinear forms, then its multi-periodic wave solutions can be directly obtained 
by solving system (10.5). 

Theorem 12. Let C(£,e',/u) and H{Dj.,Dt) be given in Theorem 10, and 
make a choice such that e'j — Sj — ±1/2, j — 1, - ■ ■ , N. Then 
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(i) If H{Dx, Df) is an symmetric operator, i. e. 

H{-D,,-Dt) = H{D,,Dt), 

then C(£, e', /u) vanishes automatically for the case when X^^^ fJ'j is an odd num- 
ber, namely 

N 

C(e, e', /Lt)|^ = 0, for ^/ij = l, mod 2. 



(ii) If H{Dx, Dt) is a skew-symmetric operator, i.e. 

H{-D,,-Dt)^-H{D^,Dt), 



then C(e,e',/u) vanishes automatically for the case when Yl!j=il-''j even 
number, namely 

N 

C(£, e', = 0, for ^A'-j = 0, mod 2. 

Proposition 12. Let — Sj = ±1/2, j = 1, • • • , N. Assume H{Dx, Dt) is 
a linear combination of even and odd functions 

H{D,, A) = H^{D,, A) + i^2(A, A), 

where ifi is even and H2 is odd. In addition, C{£,£',ijl) corresponding (10.8) is 
given by 

C(e, e', = Ci(e, e', C2(e, e', m), 

where 

Ci(e, e', /Lt) = ^ exp [-27r(T(n - /Lt/2), n - /Lt/2) - 27ri(n - /Lt/2, e' - e 

C2(e, e', /Lt) = ^ H2{M) exp [-27r(T(n - /Lt/2), n - /Lt/2) - 27ri(n - /Lt/2, e' - e 
Then 

iV 

C(e, e', /Lt) = C2(e, e', /Lt) for ^/Xj = l, mod 2, 

i=i 

C(e, e', /Lt) = Ci(e, e', /Lt), for ^/Xj = 0, mod 2. 

The theorem 2 and corollary 1 are very useful to deal with coupled super- 
Hirota's bilinear equations, which will be seen in the following section 10. 
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By introducing differential operators 



dx = Qii% + Qi2% H h QiAT^^jv = a • V, 

9* = + /52% + • • • + pNd^, = /3 • V, 



then we have 



d%^{C r) = (a • V)^(/3 • V)^^(C, r), ^, Z = 0, 1, • • • . 



10.3. One-periodic waves and asymptotic analysis 

Let us first construct one-periodic wave solutions of the NKdV equation (1.1) 
by using bilinear Backlund transformation (7.6). As a simple case of the theta 
function (10.1) with A?" = 1, s = 0, we choose F and G as follows 



where C = Cix + /St + 5 is the phase variable, and r > is a positive parameter. 

By Theorem 6, the operator Hi — D^ — X in bilinear equation (7.6) is symmet- 
ric, and its corresponding constraint equation in the formula (10.5) automatically 
vanishes for = 1. Meanwhile, H2 — DfDl — 2wqDx -\- (4i>o -|- ?>\)Dt are skew- 
symmetric, and its corresponding constraint equation automatically vanishes for 
= 0. Therefore, the Riemann theta function (10.6) is a solution of the bihnear 
equation (7.6), provided the following equations 



^{[47ri(n - ii/2)fa^ - A} exp(-27rr(n - ii/2f + m^n - I^/2)%=q = 0, 




G = ^{C, 1/2, 0|t) = J2 exp(27rm(C + 1/2) 



Trn^r) 



(10.6) 




^{[47ri(n - fi/2)fa'^P + 8m{n - fi/2)awo + 47ri(n - /i/2)(4uo + 3A)/5} 



X exp(-27rT(n - ///2)^ + TTi{n - ///2)) |^=i = 0. 



(10.7) 
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hold. 
Let 

MCp) = m, 1/4, -l/2|2r) = ^p(2n-i)^ exp[4m(n - 1/2)(C + 1/4)], 
i?2(C,p) = ^(2C, l/4,0|2r) = 5^p^"' cxp[4^7rn(C + 1/4)], 

then, the equation (10.7) can be written as a hnear system about (3 and A 

^'^a^ - i?2A = 0, 

(10.8) 



where the derivative value of 'j?j(C, p) at C = is denoted by simple notations 

dC 

It is not hard to see that the system (10.8) admits the following solution for 
the NKdV equation (1.1) 

^2 ' <'^?2a2 + 4^?i^?2^;o + 3^?;^?'2'a'■ ' 

So, we obtain the following one-periodic wave solution 

V = Vo + 2dl]n^{C,Q,Q\T), W = Wo + 2d^dtln^{C,Q,Q\T), (10.10) 

where — ax + ^t + S and parameter /3 is given by (10.9), while other parameters 
a,T,vo,wo are arbitrary. Among the four parameters, the two ones a and r 
completely dominate a one-periodic wave. 

In summary, one-periodic wave (10.10) is one- dimensional and has two fun- 
damental periods 1 and ir in phase variable ( (see Figure 3). 

In the following theorem, we will see that the one-periodic wave solution 
(10.10) can be broken into soliton solution (6.20) under a long time limit and 
their relation can be established as follows. 

Theorem 13. In the one-periodic wave solution (10.6), the parameter /3 is 
given by (10.9), other parameters are chosen as 

k 7 + vrr , , 
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(a) (b) 




(c) (d) 




-2-10 1 2 



Figure 3. One-periodic wave for the NKdV equation (1.1) with parameters: 
a = 0.6, T = 2,vo = 0.5, Wo — 1- (a) and (b) show that every one-periodic 
wave is periodic in both x and y directions, (c) Perspective view of the wave, 
(d) Overhead view of the wave, with contour plot shown. The bright hexagons 
are crests and the dark hexagons are troughs. 

where ki and 7 are the same as those in (6.20). Then under a small amplitude 
limit, one-periodic wave solution (10.10) can be broken into the single soliton 
solutions (6.20), that is, 

V — >v, W — >w, as p^O. (10.12) 

In particular, in the case of Vq = 0, Wq = 1, the one-periodic solution (10.5) 
tends to the kink-type soliton solution (5.2), that is, 

V — W — > w\ as p ^ 0. (10.13) 

Proof. Here we use the system (10.8) to analyze asymptotic properties of the 
one-periodic solution (10.10). Let us expand the coefficients of the system (10.8) 
as follows 

^[ = -47rp + 127rp9 + ■ ■ ■ , ^'1' = Wn^p + 4327rV + • • • , 

(10.14) 

= 1 + 2p^ + ■ ■ ■ , ^'^ = 327r2p' + ■ ■ ■ , 
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Suppose that the solution of the system (10.8) has the following form 
A = Ao + Aip + \2P^ H = Ao + o(p), 

(10.15) 

/5 = /5o + + (52P^ + ■ ■ ■ = /3o + o{p). 
Substituting the expansions (10.14) and (10.15) into the system (10.8) and 
letting p ^ 0, we immediately obtain the following relation 

— dWn 

— zTT^a^ + 2vq 
Combining (10.11) and (10.16) leads to 
A ^ 0, 

or equivalently rewritten as 

C = 27riC — TTT — kx + 27rif3t + 7 

2kwo (10.17) 
— ^ ~ -n : — ^ + 7 = as p ^ 0. 

It remains to verify that the one-periodic wave (10.11) has the same form as 
the one-soliton solution (6.20) under the limit p — > 0. Let us expand the function 
F in the following form 

F = 1 + p^e^"""^ + e-2"^) + p^(e^^'^ + 6"^"*^) + • • • . 

It follows from (10.11) and (10.17) 

F = 1 + + p\e-^ + e^^) + p'^e-''^ + e^^) + ■■■ 

(10.18) 

— >l + e'^ — >l + e^, as p 0. 
So, combining (10.11) and (10.18) yields 
V — y vq + 2dxx ln(l + e^), 
w — )■ Wo + 2dtdx ln(l + e^), as p — >■ 0. 
Thus, we conclude that the one-periodic solution (10.10) may go to a bell-type 
soliton solutions (6.20) as the amplitude p — > 0. □ 

10.4. Two-periodic waves and asymptotic properties 
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Let us now consider two-periodic wave solutions to the NKdV equation (1.1). 
For the case of = 2, s = 0, e = 1/2 = (1/2,1/2) in the Riemann theta 
function (10.1), we choose F and G as follows 



F = '(?(C, 0, 0|t) = exp{27ri(C, n) - 7r(Tn, n)} 



G = ^?(C, 1/2, 0|t) = J] exp{27ri(C + 1/2, n) - n{Tn, n)} 

= (-1)"^+"^ exp{27ri(C, n) - 7r(Tn, n)} 
Tiez2 



;io.i9) 



where n = (^1,77.2) G ^^ C = ((1,(2) e Q = ajx + f3jt + 6j, j = 1,2, 
and cx = (01,0:2), /3 = (/9i,/92) G C^. The matrix r is a positive definite and 
real- valued symmetric 2x2 matrix that is, 

T = (%)2x2, T12 = T21, Til > 0, T22 > 0, T11T22 - T12 > 0. 

According to Theorem 5, constraint equations associated with Hi — Dl — X 

and H2 = DtDl — 2wqDx + {^vq + 'i\)Dt automatically vanish for {^1,^2) = 
(0, 1), (1, 0) and for (/ii, /U2) = (0, 0), (1, 1), respectively. Hence, making the theta 
functions (10.19) satisfy the bilinear equation (7.6) gives the following constraint 
equations 

^ [-167r^(n - /u/2, a)^ - A] exp{-27r(T(n - /u/2), n - /u/2) 

2 

+ m ^{nj - /Xj/2)}|^=(^i,^2) = 0, for (//i, = (0, 0), (1, 1) = 0, 

i=i 

^ [-647r^i(n - /Lt/2, a)^(n - /u/2, /3) + 87ri(n - /u/2, a)wo + 47ri(n - /u/2, /3)(4t;o + 3A)] 

2 

X exp{-27r(T(n - /u/2), n - /x/2) + ni ^{nj - /ii/2)}|;,=(/,,,^2) = 0, 
for (//i,//2) = (0,l), (1,0). 

(10.20) 
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Let 

Pki = 6"'^""'='/^ k,l = 1,2, p= (pii, pi2, P22) 
^,(C,p)=^(2C,l/4,-V2|2r) 

ni,n2€Z k,l=l 
Sj ^ {Sj,l,Sj,2), J = 1,2, Si = (0,1), S2 = (l,0), S3 = (0,0), S4 = (l,l) 

then the system (10.20) can be rewritten as a hnear system 

(a • Wf'dj - X'dj = 0, j = 3, 4, (10.21) 

(/3 ■ V)(a ■ V)% + 2wo{a ■ + {Avq + 3A)(/3 ■ = 0, j = 1, 2, 

(10.22) 

where "dj represent the derivative values of functions "^jiC, p) at Ci = C2 = 0. 
The system (10.22) admits a unique solution 

_ \ dif,g) 1 ' f2wo{a ■ V)7?i\ . 

W-[a(Ci,C2)J [2wo{cx-v)^2j ^ ■ ^ 

where 3/'^^^ is the Wronskinan matrix given by 



9(Ci,C2) v^ci^ 

/ = [(a • V)' + 4t;o + 3A]^?i, g = [(a • V)' + Av^ + 3A]i?2. 
With the help of the above (^1, ^2), we are able to get a two-periodic wave solution 
to the NKdV equation (1.1) 

V = vo + dl\n^{CQ,Q\T), W = wo + d,dt^{CQ,Q\T), (10.24) 

where ai,a2, T12, 5i and 82 are arbitrary parameters, while other parameters ^2 
and Til, ^22 are given by (10.23) and (10.21), respectively. 

In summary, the two-periodic wave (10.24) is a direct generalization of two 
one-periodic waves. Its surface pattern is two-dimensional with two phase vari- 
ables Ci and C^2- The two-periodic wave (10.24) has 4 fundamental periods {ei, 62} 
and {iTi,iT2} in (Ci,C2), and is spatially periodic in two directions Ci;C2- Its real 
part is not periodic in 9i direction, while its imaginary part and modulus are all 
periodic in both x and t directions. 
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Figure 4. Two -periodic wave for the NKdV equation (1.1). (a) and (b) 
show that every one-periodic wave is periodic in both x- and y-directions. (c) 
Perspective view of the wave, (d) Overhead view of the wave, with contour 
plot shown. The bright hexagons are crests and the dark hexagons are troughs. 



Finally, we study the asymptotic properties of the two-periodic solution (10.24). 
In a similar way to Theorem 5, we figure out the relation between the two-periodic 
solution (10.24) and the two-soliton solution (6.21) as follows. 

Theorem 14. Assume that {Pi, P2) is a solution of the system (10.22), and 
in the two-periodic wave solution (10.24), parameters aj,Sj,Ti2 are chosen as 



I3 + TTT- 



33 



A 



ri2 



12 



J = 1,2, 



(10.25) 



27rz 27r 

where kj,'^j,j = 1,2 and A12 are those given in (6.21). Then, we have the 
following asymptotic relations 



A - 
F 



0, 



HT. 



33 



3 



1,2, 



[10.26) 



1 + + e''^ + e''l+''2+^l^ as pu,p22 ^ 0. 

So, the two-periodic wave solution (10.24) just tends to the two-soliton solution 
(6.21) under a limit condition 



V 



V, W 



as pii,p22 0. 
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Proof. Using (10.20), we may expand the function F in the following exphcit 
form 

F =1 + (e^^^^i + e-2^*^i)e-^^ii + (e^^^^^ ^ g-2«C2)g-7rT22 

_^^g27ri(Ci+C2) _^ g-27ri(Ci+C2)^g-7r(Tii+2ri2+T22) _^ 

Furthermore, adopting (10.25) and making a transformation we infer that 

F = 1 + e^i + e^^ + e^i+^^-'^^^i^ + pt^e'^' + p^^e'^^ + pt^pi^e'^'-^'-^''^'' + ■■■ 
^ 1 + e^i + + e^^+^2+^l^ as pii,p22->0, 

where Q — ajx + j3jt + j = 1, 2, and $j — 27riPj,j — 1,2. 
Now, we need to prove 

^2^4^^ ^ Cj, J = 1,2, as pii,p22 ^ 0. (10.27) 

As in the case of A?" = 1, the solution of the system (10.23) has the following 
form 

Pi = A,0 + + ^2,2P22 + 0(pii, P22), 

/32 = /32fi + /32,iPn + /32,2P22 + o(pii, P22), (10.28) 

A = Ao + Aipii + A2P22 + o(pii, P22)- 
Expanding functions "dj, j = 1,2,3,4 in equations (10.21) and (10.22) with 
substitution of assumption (10.28), and letting pn, P22 — > , we will obtain 

Ao = 0, 

167ri(-7r^ai + t;o)/5i,o - Sniwoai = 0, (10.29) 

167ri(— TT^tta + '^o)/32,o — STiiwoa2 = 0. 
Using (10.28) and (10.29), we conclude that 

A = o(pii,p22) — > 0, 

—2kjWo —2kjWo 
= + o(pn, P22) ^ as pn, P22 ^ 0, 

and therefore we have (10.26). So, the two-periodic wave solution (10.24) tends 
to the two-soliton solution (6.21) as pn, P22 — >■ 0. □ 

10.5. Multi-periodic wave solutions 
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The system (10.5) indicates that constructing multi-periodic wave solutions 
depends on the solvability of the system (10.5). Obviously, the number of con- 
straint equations of the type (10.5) is 2-^"^ -|- 1. On the other hand, we have 
|iV(iV + 1) -|- 3A^ -|- 3 parameters Tjj, Tjj, a^, cuj. A, Mq, Vo- Among them, 2A'" pa- 
rameters Tii, ai may be the given parameters relate to the amplitudes and wave 
numbers of A^-periodic waves. Therefore, the number of the unknown parameters 
is \N[N -|- 1) -|- a*" -|- 3 while }^N{N + 1) parameters r^j, implicitly appearing in 
the series form, can not to be solved explicitly in general. So, the number of 
the explicit unknown parameters is only A^ -|- 3, and the number of equations is 
larger than the unknown parameters in the case of A^ > 4. This fact means that 
if equation (10.5) is satisfied, then we have at least A^-periodic wave solutions 
(A" < 4). In this paper, we only consider one- and two-periodic wave solutions of 
the NKdV equation (1.1). There are still certain computation difficulties in the 
calculation for the case of A^ > 2, which will be studied in the future. 
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